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Abstract. We study the energy flow of dissipative dynamics on infinite lat- 
tices, allowing the total energy to be infinite and considering formally gradient 
dynamics. We show that in spatial dimensions 1,2, the flow is for almost all 
times arbitrarily close to the set of equilibria, and in dimensions > 3, the size 
of the set with non-equilibrium dynamics for a positive density of times is 
two dimensions less than the space dimension. The theory applies to first and 
second order dynamics of elastic chains in a periodic or polynomial potential, 
chains with interactions beyond the nearest neighbour, deterministic dynamics 
of spin glasses, discrete complex Ginzburg-Landau equation, and others. We 
in particular apply the theory to show existence of coarsening dynamics for a 
class of generalized Frenkel-Kontorova models in bistable potential. 



1. Introduction 

The state space of lattice dynamical systems we study consists of functions u : 
— ^ R*^. Here N denotes the dimension of a lattice, and each lattice point has 
M degrees of freedom. We develop here general results on asymptotic behavior 
of lattice systems which can be understood as a thermodynamic limit of a family 
of finite-dimensional, gradient dynamical systems. The motivating example is the 
Frenkel-Kontorova (FK) model, with the total energy given formally by 

oo 

E{u) = - + V{u,), 

■i— — oo 

where is a smooth periodic potential and u = (ui)igz G Its gradient dynamics 
is given by dui/dt — d / duiE{u) (see Section 3 for details). The FK model has been 
studied as a reasonable physical model (e.g. for Josephson junction arrays) with 
rich features present in more dimensional and more complex systems (see [2], [3], 
[TU] . [T5] . [18] , pp] and references therein). The FK model is an example 

of three paradigms studied in this paper: this is an example of an extended, a 
dissipative, and a system on a lattice. 

The extended systems property in the setting of partial differential equations 
means that a particular state of the system is defined on a unbounded set (typi- 
cally the entire M^), with no required convergence to at infinity. Their physical 
importance is that such a large state space is required to contain solutions such 
as travelling waves, kinks and multikinks etc. In the setting of lattice systems, we 
consider systems of infinitely many ordinary differential equations, and there is no 
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straightforward reduction to a finite dimensional system (e.g. because of periodicity 
of solutions). 

We focus here on the dissipative structure of extended systems. For (unforced) 
systems for which the total energy is finite, the dissipative structure typically im- 
plies existence of a Lyapunov function - typically the energy functional itself - which 
is strictly decreasing along any non-stationary solution, and bounded from below. 
By the LaSalle principle, this implies a complete asymptotic description of the dy- 
namics: the dynamics asymptotically converges to the set of equilibria. Also strict 
local minima of the Lyapunov function are Lyapunov stable. For extended dissi- 
pative systems, the techniques do not apply as the energy functional is divergent. 
As we will see, the dynamics resembles gradient systems, but can and typically is 
topologically different. 

Th. Gallay and the author have explored the extended dissipative structure 
for partial differential equations. In [T3], [11], we have shown that a large class 
of extended dissipative PDEs, including examples such as the reaction-diffusion 
equation and the Navier-Stokes equation, have certain universal features. In short, 
for spatial dimensions 1,2, the dynamics is for almost all times arbitrarily close to 
the set of equilibria. In dimensions iV > 3, the persistent non-equilibrium dynamics 
can be at most of the dimension R^^^. The dynamics, however, topologically 
typically differs from gradient systems, as in dimensions iV > 1 w-limit sets can 
contain non-equilibrium points, and in dimensions TV > 3 periodic orbits and more 
complex recurrent behavior can exist. 

In this paper we explore this in the context of lattice systems. We show that 
the same essential conclusions hold, but with modified proofs and the key energy 
flux and dissipation estimates. We also show that the conclusions apply to several 
examples which are difficult to model as PDEs, such as lattices with interactions 
beyond the nearest neighbour, and lattices with randomly chosen interactions for 
each pair of lattice points, such as spin glasses. In [21], we already proved some of 
the qualitative features for lattices in dimension iV = 1. We here also strengthen 
the results from [21j even in dimension iV = 1, as we provide quantitative bounds 
related to all claims, whereas the results in |21| are merely of qualitative nature. 

Our approach is axiomatic: we specify universal features of the energy flow of 
dissipative lattice dynamics in Section 3, and then deduce general theorems. As 
in |14l . we assume a solution of a system of equations is a semiflow ip on a state 
space X, and then assign to each u € X its energy, dissipation and flux functions 
e,d,f. Whereas the Lyapunov function associates to each u€ X a single number, 
in our approach we associate to each u e X its energy, flux and dissipation at each 
point a £ of the lattice. We call such systems the Lattice Extended Dissipative 
Systems, or shortly Lattice EDS. The key features are the energy balance inequality, 
and the property that flux generates dissipation. We express it by using the notation 
analogous to calculus on M^, adapted to the discrete space Z^. In Section 4 we 
then show that many examples satisfy our assumptions, including elastic chains in 
any dimension and with various potentials, but also systems already mentioned, 
and semi-discretizations of PDEs such as the discrete complex Ginzburg-Landau 
equation. 

In Section 6 we show that for dimensions N — 1,2, we can find R large enough 
(as a function of T), so that the total energy is not increasing on the cube of radius 
R (this is not necessarily true if iV > 3). The proof of this encodes the key ideas 
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and reveals key constraints to the energy flux for Lattice EDS. In Sections 7 and 8 
we develop universal bounds on the energy flux and dissipation, which we then use 
in Section 9 to describe the asymptotics of Lattice EDS. We also use the flux and 
dissipation bounds to develop general bounds to relaxation times. 

In Section 11 we consider an application of the results to coarsening in a bistable 
potential. This has been studied for continuos-space systems by Eckmann and 
Rougemont in dimensions = 1 and iV = 2 ([7], [H]). In particular, they have 
shown that for the real Ginzburg-Landau equation, topologically a continuous ana- 
logue to gradient FK model dynamics without forcing, one can construct orbits 
such that the solution u{x) at each point x of the physical space "jumps" infinitely 
many times between two stable equilibria. We call this phenomenon coarsening 
dynamics. Here we study a class of FK models in dimensions N = 1,2. We give 
sufficient conditions on a probability measure /i on the state space, so that /i-a.e. 
initial condition u : — > R has coarsening dynamics. This means that for each 
lattice point a G , u{a) jumps infinitely many times between two stable equi- 
libria of the potential as f — oo. We show this for both the first degree, gradient 
dynamics without forcing (analogously to Eckmann and Rougemont), but also for 
the second degree dynamics with sufficiently strong damping. 



2. The calculus on lattices 

In this section we recall the notation and elementary results of calculus on lat- 
tices. We use the analogy between the standard results such as partial integration 
and the Stokes theorem, with minimal technical difficulties due to discreteness of 
the space. We denote by u,v,w the functions ^- M or ^ R^^ The letter N 
will always denote the dimension of the lattice {d being reserved for dissipation). 
Greek letters a, /3 will denote the elements of Z^. Let Sj G be the "basis" for 
j = 1, N: {ej)i = for i j, 1 for i = j. The translation Tj is defined naturally 
for u : Z^ R or w : Z^ R*^ with Tju{a) = u{a + 6^). 

The partial derivatives and differentials for u : R are now: 

dju{a) = u{a) — u{a — Ej), 

d*u{a) — u{a + ej) — u{a) — Tjdju{a), 

'S/u{a)j — dju{a), 

V*M(a)j = d*u{a), 

where 'Vu{a), 'V*u{a) : Z^ — R^. For = 1, we will also use the notation (5* . 
Given v : Z^ — >■ R^, the divergence is 

N 

divu(Q!) — djv{a)j, 
i=i 

N 

div* (a) = 
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We can also naturally define the Laplacian Au : — > M for u : — > M by any 
of the following equivalent relations: 

(2.1) Au = div V*u = div* Vm 

N N 

(2.2) = ^a,5> = ^9;9,u, 

3=1 j=i 

or equivalently Au{a) = '^jLi [u{a + Sj) + u{a — Sj) — 2u{a)]. The proof of the 
following partial integration relations is straightforward and left to the reader. 

Lemma 1. Let u, w : R. Then 

(2.3) d*{uw) = d*u-TjV + u-d*v 

(2.4) = dju ■ V + TjU ■ djv, 

(2.5) div*(u-Vw) ^ uAw + V*u-V*v, 

(2.6) div(u-V*w) = uAw + Vu-Vv. 

Our key tool will be the analogue of the Stokes theorem. For square lattices, the 
natural setting is to consider the cubes C{r),C*{r) of the radius r gN, and their 
boundaries dC{r),dC*{r), defined as follows: 

C{r) = {aGZ'^, - r <aj <r-l, j = 1,...,N}, 
C*{r) = {aGE^, -r + l<aj <r, j = 1,...,N}, 

dC{r) = dC*{r) = {a G Z^, \aj\ < r for j = 1, ...,7V, but for some k, \ak\ = r}. 
The normal vectors n : dC{r) — >■ M^, n* : dC*{r) — > are now 

X^iLi ^ji ) c^ii ) ctjfc = '^1 all other — r + \<aj<r— 1 
n*{a) = ^ — ej. au = —r, all other — r + l<aj<r 

otherwise, 

£fe ak = r, all other — r<aj<r — 1 

n{a) = ^ -J2i=i£ji, aji, = -r, all other - r + 1 < aj < r - 1 
otherwise. 



Note that |n*(a)|, |n(a)| < vN, with the values higher than 1 occurring at 
"corners" (at which also n*{a),n{a) differ). 

We push further the analogy with the continuous case, hopefully to ease reading 
of the arguments to follow, by often using the integral symbol instead of the finite 
sum over a bounded set A c Z^: 



u{a)da := ^ u{a). 

The discrete space analogue of the Stokes theorem now reads: 
Lemma 2. Let v:Z^ andr>l. Then 



(2.7) / divv{a)da = / v{a) ■ n*{a)da, 

Jc'(r) JdC*{r) 

(2.8) / div* v{a)da = I v{q) ■ n{a)da. 

Jc(r) JdC(r) 
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Proof. First note that by definition 

a6zC*(r) QGC*(r),aj— r aGC* (r) ,aj — — r+1 

(all the other terms cancel out). We now sum over j = 1, r, recall the definition 
of n* and get (PTT)) : is analogous. □ 



3. The key properties 

In this section we specify key properties of the energy flow of dissipative systems 
on lattices. To motivate the definitions to follow, we start with perhaps the simplest 
example: the over-damped dynamics of the one-dimensional FK model, i.e. a 
system of elastically connected balls in a periodic potential. A state of the model 
is given by a function u : Z R, and its total energy is a formal sum (typically 
divergent): 

(3.1) E{u) = - - 1) - + ^(«(")) 

where 6 G M is a constant (the characteristic length) , and : K — > M is a smooth 
periodic potential, V{x + 1) = V{x). The gradient dynamics is formally expressed 
as 9(U = — Vi?(u), or more precisely with the equations 

(3.2) dtu{a) = -2u{a)+u{a-l) + u{a + l)-V'{u{a)) 

= l\u{a) -V'{u{a)). 

It is well known that the equation p.2p defines a continuos semiflow ip on the 
space Xk = {u : Z — >■ R, |u(q;) — u{a — 1)| < K for all a}, where K is any positive 
integer. In particular, local existence and continuity follow from the standard results 
on solutions of ordinary differential equations on Banach spaces, and Xk is invariant 
by the order-preserving property of (13. 2[) and the existence of stationary solutions 
with the spacing K . We consider two topologies on Xk'- the first one is induced by 
the 1 1. 1 loo norm, and the second one is the induced product topology (as a subset 
of R^). The semifiow ip is continuous in both topologies. 

As was noted in the introduction, the dynamics (j3.2p is not gradient in the strict 
sense (and the LaSalle principle does not hold), as the "Lyapunov" function (13. ip 
is typically not finite. 

We can associate to each -u : Z — > R its local energy e„(a) as 

e„(a) = i(u(a) - u{a - l)f + V{u{a)), 

and then the total energy in the interval [— r + l,r] = C*(r), denoted by Eu{r) is 
equal to (using also the notation from the previous section) 

Eu{r)= f eu{a)da= V \]-{u{a) ~ u{a ~ l)f + V(u{a))\ . 
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We now consider the time evolution of the energy Eu{r), and obtain the energy 
balance equation: 

r 

dtEu(t){r) = {-{Mo)-u{a-l)-u{a + l))dtu{a) + V'{u{a))dtu{a)} 

a— — r+1 

~(u(r + 1) — u(r))dtu{r) + (u{—r + 1) — u(—r))dtu{—r) 

r 

= - ^ {dtu{a) f ~ d*^u{r)dtu{r) + d*au{-r)dtu{~r). 

a— — r+l 

We can now introduce the local energy dissipation d„ and flux associated to 
M : Z ^- R with 

du{a.) = {dtu{a)Y, 

fu{a) = -dtu{a){u{a + I) ~ u{a)) ^ dtu{a)dl^u{a). 

The energy balance equation can now be written in its integral and local form, also 
related by the Stokes theorem (Lemma [2]) 

(3.3) 9* (r) = I d^(t){oi)da+ I fu(t){Q) ■ n* {a)da, 

JC*(r) JsC-{r) 

(3.4) dte = -d + div/. 

The energy balance equation p.4p is the first key component of the theory to 
follow. The second one is the observation that there is no energy transport (or 
energy flux) without dissipation. For the equation p.2p we can without loss of 
generality assume V > (otherwise we add a constant). With the choice of e,d,f 
as above, we have. 

f < 2||e||ood. 

Motivated by this, we introduce an abstract definition of an extended dissipative 
system on a lattice, associated to a continuous semifiow rather than to a specific 
set of equations. If is a semiflow on a state space X, we will associate to each 
u € X the functions describing energy, dissipation and flux at each point of a lattice 
Z^. Thus we have a family of mappings e, d, /i, /jv : X — )• /oo(^^), where 
/ = (/i, /at) is the flux vector. Here Zoo(^^) is as usual the space of bounded 
functions u : — >■ R with the sup- norm ||.||oo- A point u ^ X is stationary for 
the semiflow ip, if for each t > 0, (p{u, t) = u. 

We now list the required properties of the functions e, d, f, as we have seen 
satisfied for the equation p.2|) . 

(Al): e,d> 0; 

(A2): For each a S Z^, the functions X i-7> e{a),d{a), fi{a), fN{ct) are 
continuos as functions X R; also t ||e„(t)||oo, |Mti(t)||oo, ||/«(t)||oo are 
continuous as functions R+ R+; 

(A3): If d = then u G X is stationary; 

(A4): The energy balance equation holds (dt is derivative in terms of distri- 
butions): 

(3.5) dte^-d + divf] 

(A5): There is a non-decreasing function b : R+ R+ such that 

(3.6) <6(||e||oo)-d. 
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Definition 1. Assume X is a metrizable space and ip a continuous semiflow 
on X. We say that {X,(p) is a N- dimensional extended dissipative system on 
a lattice (Lattice EDS), if there exist functions e, d, /i, /jv : X — >■ ^oo(Z^); 
/= satisfying (Al)-(A5). 

In the dimension iV = 1, we will sometimes write (|3.5|) as 

dte = -d + daf- 

Given a point u G X, wc will often denote its semiorbit by u{t), and the family 
of functions e,d,f with e{a^t),d{a,t)^ f{a^t) (evolution of the energy, flux and 
dissipation along a semiorbit). Note that then (A2) and continuity of the semiflow 
imply that e, d, /i, /at : If^ x [0, oo) — >■ R are continuous in the second variable. 
When not clear to which u ^ X we associate e,d,f, we write it as e„,(i„,/„, 
otherwise the subscript u is omitted. 

In the next section, we give a number of examples justifying the definition of 
Lattice EDS. 

We briefly comment perhaps the most technical assumption (A2). It may seem 
that this can be replaced with a simpler requirement that e,d, f are continuous as 
functions X looC^^)- This will however not typically be true in our setting, as 
we will typically consider the product topology on X (i.e. topology of pointwise 
convergence at lattice points). Thus we require only that e,d,f are continuous 
locally (i.e. at one lattice point), and uniformly continuous along a semiorbit, 
which will hold in our examples. 

In most cases, we will also require that the energy e is uniformly bounded on X , 
so we introduce the following notion: 

Definition 2. We say that a Lattice EDS is bounded, if sup ||eti||oo < +oo. 

For bounded Lattice EDS, we often denote by /? = sup^gj,;- ^(||eu||oo), and then 
(|3.6|) becomes 

(3.7) /2 < 13d. 

This requirement means that the local energy remains uniformly bounded, while 
the total energy may still be (and typically is) infinite. This is typically physically 
a reasonable requirement for dissipative systems, and in a number of cases can 
be demonstrated by either the ordering property of solutions (i.e. a maximum 
principle), or by appropriate local energy estimates. We will discuss it further 
related to particular examples. 

4. Examples of extended dissipative systems on lattices 

4.1. The standard damped Frenkel Kontorova model. The set of states of 
a N-dimensional FK model is the set of functions / : — > R^. The dynamics 
is governed by the elastic force between the particles moving in a smooth periodic 
potential V : — >■ R, V{x + a) — V{x) for all a G Z^. The damped equations 
can be thus written (by again relying on the notation from Section 2): 

(4.1) A • dttujia) + dtujia) = Au,{a) - — l/(u(a)), 

where u : -> R^, u{a) = (Mi(a), U2(a), ujv(a)), and A > is a constant. 
Again without loss of generality we assume V > 0. We define the set of functions 
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/ : ^ of bounded width with 

(4.2) Y ^{u:Z^ ^ K", ||Vw||oo < oo}, 

where Vu(a) is a x iV matrix, {Wu{a))i,j = diUj{a), its norm |Vu(a)p = 

a-iid l|V"lloo = sup^gzjv |Vm(q!)|. 

In the gradient case A = 0, the state space will he X ~ Y, and in damped the 
case A > 0, it will be X = F x /oo(Z^), {u,dtu) e X. It is well known that the 
equation (|4.ip generates a continuos semiflow on X, The considered topology is the 
induced product topology on Y (as a subset of (M^)^ ) and in the case A > the 
product topology on X. In particular, the global existence of solutions and the 
invariance of X follows from the fact that ||Vu||oo can grow at most exponentially, 
thus it can not diverge at some finite time. 

We now show that this is a N-dimensional Lattice EDS. We define the energy, 
dissipation and flux associated to u G X with 

e(a) = ^idMc^)f + \\Vuia)\^ + Viuia)), 
d{a) = {dtu{a)f, 

N 

where [dtu)"^ = J^jLii'^tUj)'^- We easily see that in both cases A = 0, A > we 
have e,d, f € ^oo(Z^), and that the axioms (Al), (A2), (A3) hold trivially. We now 
have 

N N ^ d 

dte{a) = ^ \dtUj ■ duUj + ^ Vuj(q;) • dtVuj{a) + ^ -^V{u{a))dtUj{a). 

As by the partial integration (j2.6|) . the second sum on the right-hand side is equal 
to X]j=i div(9tUj(a) • V*Mj(a))— dtUj(a)l\uj(pL), we get the energy balance 

equation 

N N 

dte{a) = -^(atU,(a))2 + ^div(9tu,(a)-V*Mj(a)) 

= — (i(a) + div /(a), 

thus (A4) holds (note that e is differentiable for all t > Q). We can also easily check 
(A5): 

(N \^n/n 

N N 

< Y.(9tuj{a)rj2{d:u,{a)r < d(a)|V*^.(a)p 

1 2 

< d{a) sup |V*M(a)p < d(a)— sup |Vu(a)p < — d(a)||e| 

aGZ« A* c<ez« A 

We conclude that ()4.ip generates a Lattice EDS on A. 
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One can consider more specific state spaces. Let us first discuss the case A = 0. 
For a given constant K > 0, let Yk be the set of all states whose semiorbits have 
uniformly bounded width 

Yk^{u:Z^ ^ R^, ||Vu(t)||oo < K for ah t > 0}. 

Then by definition, Yk is invariant, and (|4.ip generates a Lattice EDS on Yk too. 
We already noted that ior N — 1, all u ^ Y are in Yk for a sufficiently large 
integer K > 0. Also in any dimension, all periodic orbits of ()4.ip are in Yk for K 
large enough. Also e on Yk is bounded, thus the Lattice EDS is bounded. If we 
introduce the equivalence relation U'^viiu — v = a€z Z^, then one can show 
that the dynamics (14.11) is well defined on the quotient space Yk/ ^, and induces 
a Lattice EDS on Yk/ ^ (the definitions of e,d,f are invariant with respect to 
~). By TychonofF theorem, Yk / ~ is compact in the induced product topology. 
Thus the study of (|4.ip reduces to a study of a semiflow on a compact space with 
a Lattice EDS structure. 

An analogous construction is possible also in the case A > 0, by also uniformly 
bounding the first derivatives dtu. 

4.2. Generalized Frenkel-Kontorova model. We now discuss chains with more 
general nearest-neighbour interactions. We consider one-dimensional chains in a 
larger-dimensional space, i.e. the set of states comprises of functions u : Z — > M*^. 
We assume that the nearest-neighbour interaction is defined by a smooth potential 
L : K^^ X M*^ M satisfying 

(LI): L{x, y) = L{x + a, y + a) for all a G Z^, 
(L2): YmY\^_y\^oo L{x,y) = +0O. 

Clearly the standard ID Frenkel-Kontorova model is a special case of this. The 
gradient equations of motion are 

(4.3) dtu{a) = —L2{u{a — l),u{a)) — Li{u{a), u{a -\- 1)), 

where Li, L2 denote partial derivatives with respect to the first and second coordi- 
nate. We consider the gradient dynamics only, though the damped one is also an 
example of a Lattice EDS (the choices of e, d, / and calculation are analogous to 
the standard case). 

We note that this example includes a particularly important family of potentials 
L: generating functions of M-dimensional Tonelli Lagrangian maps. These are 
discrete analogues of Tonelli Lagrangian flows on a torus (see e.g. [5]), and can 
also be understood as time-one maps of Lagrangian flows. In addition to (LI), 
(L2), Tonelli Lagrangians satisfy the following: L12 is positive definite, and (L2) is 
strengthened to be lim|2._j^|^oo L[x, y)/\x — y\ = +00. Then the stationary points 
of (|4.3p correspond by variational principle to orbits of the associated Lagrangian 
map, and the study of (j4.3p is the study of formally gradient dynamics of the action 
functional. 

Without loss of generality L > (by (LI), (L2), L is bounded from below; so if 
not we add a constant). We define e, d, / associated to a state u : Z — M*^ with 

e{a) = L{u{a — 1) , u{a)) 

d{a) = {dtu{a)f', 

f{a) = —Li{u{a),u{a + l))dtu{a). 
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We first check (A4): 

dte{a) = Li{u{a — l),u{a))dtu{a — 1) + L2{u(a — l),u{a))dtu{a) 

= —{dtu{a))'^ + Li{u{a — 1), u{a))dtu{a — 1) — Li{u{a), u{a + l))dtu{a) 
= -d{a) + daf{a). 

The property (A5) will follow directly from this: 

Lemma 3. Assume L : continuously differ entiahle and satisfies 

(LI), (L2). Then there exists a non- decreasing function b : K+ — > M+ such that for 
all x,y (E R*^ 

L,ix,yf <bi\Lix,y)\). 
Proof. We define the following functions 6i , 62 '■ ~> as 

bi{x) = sup Li{a,b)^, 

|6— a| <x 

hiy) = sup \b-a\. 

\via.b)\<y 

We see that &i is < 00 by periodicity (LI): the domain of the supremum can 
be reduced to a compact set. Also 62(2/) < 00 because of (L2). By definition 
61,62 are non-decreasing and Li(a,6)^ < 6i(|6 — a|), |6 — a| < 62(1^^(0,6)!). We 
set 6 = 62 o 61 , and finally note that the composition of non-decreasing functions is 
non-decreasing. □ 

Again let X be the set of states u : Z — > M.^^ of bounded width as in (14. 2p . with 
the induced product topology. By combining all of the above, we see that (|4.3p 
induces a continues semifiow on X which is a Lattice EDS. 

4.3. Lattices with interactions beyond the nearest-neighbour. For simplic- 
ity, we discuss one-dimensional lattices with the set of states u : Z — > M in the 
gradient case, though higher dimensions and the damped case can be analysed 
analogously as above. Assume C C Z is a finite set of lattice distances for which 
there is a non-zero interaction potential L'' : R x R — )■ R, 7 € C, satisfying (LI), 
(L2) and without loss of generality L'^ > 0. For simplicity assume ^ C (if there 
is a periodic site potential, we can add it to any L'^). The equations of motion are 
now 

(4.4) dtu{a) = ~Y^ {Ll{u{a - 7), u{a)) + Lj{u{a),u{a + 7)). 

-yec 

Let X be again the set of all w : Z — > M of bounded width, and let e,d, f be 
d{a) = {dt{a))\ 

/(") = -^Ll{u{a),u{a + ^))dtu{a). 

The proofs of (Al), (A2), (A3) are straightforward, and (A4) is analogous to 
previous examples. 
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Assume 6''', 7 G C are functions constructed in Lemma |3l associated to L'^ . Then 
f{a) = {dtu{a)f\Y,Ll{<a),u{o^ + ^))^ 

< |C|(9,u(a))^^6^(||e|U). 

We see that (A5) holds with b — \C\ J2^£C ■ Finahy we note that if for ah 7, 
L'^ satisfies the twist condition ^ Oi then the equation (|4.4p preserves ordering. 
We can then show that Xk — {\\dau\\ < K} is invariant for aU integer K, and (|4.4p 
generates a bounded Lattice EDF on Xk- 

4.4. Spin glasses and other random interaction models. We can also define 
a lattice in which the interactions are chosen from a finite family of interactions 
satisfying (LI), (L2), randomly for each pair of neighbouring lattice points. We give 
an example of damped dissipative dynamics of a system which is closely related to 
the Edward- Anderson spin-glass model though numerous other variations are 
possible. Assume the site potential is given by 

It is bistable, with minima at a; = ±1. Choose parameters /x, A > 0. We define two 
possible nearest-neighbour interactions 

L-{x,y) = ^(x-y-2f. 

As in spin-glasses models, the interactions either favour the lattice points being in 
the same, or the opposite site potential minima. Assume the set of lattice states 
is the set of functions m : — ^ M. Let S : 1? y. {— ei,ei, —£2, £2} {+, — } be a 
random function assigning to each neighbouring lattice pair a symmetric interaction 
(i.e. satisfying S'(a, — £j) = S{a — ej,ej)). The damped lattice dynamics is then 
given by 

(4.5) X ■ dttu{a) + dtu{a) ^ ~ ^ L^'^"''\uia),u{a + e)) + u{a) - u{a)^ . 

e— ±ei ,±£2 

Then (|4.5p defines a continuous semiflow on X = {u : — ^ K, ||m||oo < !}• 
Because of the ordering property of (|4.5p . X is invariant, and it is compact in the 
induced product topology. Analogously to previous examples, we can easily show 
it is a bounded Lattice EDS, with e, d, f chosen as follows: 

e(a) - ^(dtuf + l E L^^''^'\<c.),u{a + e)) + V(u(a)), 

e=±ei ,±£2 

d{a) = {dtu{a)f, 

/(") = E Lf''"''''\u{a),u{a + e))dtu{a). 
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4.5. The discrete complex Ginzburg-Landau equation. The discrete Ginzburg- 
Landau equation is the spatial discretization of the well-known complex Ginzburg- 
Landau equation ([I]), also studied independently, with solutions the time evolu- 
tions of functions u : l/' — > C. The equations are given with 

(4.6) ut{a) = (1 + i\)/Su + u - (1 -I- i\)\u\^u, 

where A > is a real parameter. We consider the dissipative case (in the general 
case, the two occurrences of the constant A in ()4.6p are not necessarily the same). 
Following "T^, we define the auxiliary function v{a^t) = u{a,t)e'^^^ , and then 

(4.7) vt{a) ^ {l+i\){^v + v ~\v\^v). 

One can check analogously to previous examples that (|4.7I) generates a Lattice EDS 
with 

e{a) = i|V«(a)p + i(l-|«(«)n2, 



d{a) 



\vtia)\' 
l + A^ ' 



/(a) = Re{yt{a)V*v{a)^ 



One can show, by adapting the local energy estimates to the discrete case, that it is 
also (by choosing a bounded invariant set X) a bounded Lattice EDS on a compact 
set. 



5. Energy increasing solutions of Lattice EDS 

The fundamental property of a gradient systems is that its energy functional (or 
Lyapunov function) E{u{t)) strictly decreases along any non-stationary orbit. A 
natural question to ask in the framework of extended dissipative systems on lattices 
is: are there any constraints on the set of integers i?, for which 



E{R,T) :^ / e{a,T)di 



is for a given T > not smaller than E{R,0). As already noted, we associate the 
functions e,d, f : x [0, T] to a fixed orbit u{t) of a semiflow in X. 

The key tool in the following will be the integral form of the energy balance 
equality (A4), which follows from (A4) by the Stokes theorem and integration over 

[o,r]: 

(5.1) 

f{a, t)n* [a)dadt. 

'(Ft) Jo JdC'iR) 

We introduce the notation for the total dissipation in the ball C*{R), and the 
total flux through dC*{R) up to the time T > 0: 

(5.2) D[R,T) = [ [ d{a,t)dadt, 

Jo Jc*{R) 

(5.3) F{R,T) ^ [ [ f{a,t)n*{a)dadt. 

Jo JdC'(R) 



E{R,T)^E{R,0)^~ [ [ d{a,t)dadt+ [ [ 

Jo Jc'(R) Jo JdC 
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Now (|5.1|) can be written as 

(5.4) F{R, T) = E{R, T) - E{R, 0) + D{R, T). 

We will use that to prove the following restriction on the energy increasing solutions 
in dimensions N = 1,2. We denote by Jt the following set 

JT = {Re N, E{R, T) > E{R, 0)}. 

Theorem 1. Assume {X,ip) is a Lattice EDS, and choose non- stationary u € X. 

(i) If N = 1, then Jt is finite. 

(ii) If N = 2, then X^rG/r finite. 

The claim will follow from the following Lemma, whose proof is essentially a 
discrete analogue of the separation of variables technique. 

Lemma 4. Assume J = {rp < ri < ...} is a set of positive integers, and assume 
there is a recurrent sequence satisfying the following for some A > 0.' 

(5.5) Go = e>0, 

(5.6) Gfc = Gfe_i + AGti/rf_Y, 

(5.7) Gl^i/rj^Si is bounded. 

If N = I, then J is finite. If N = 2, then X^re/ 1/^ '■^ finite. 
Proof. As Gfc > 0, dividing 

(ESI) by G^._i we get 

A _ Gfc 1 

rf_V Gu-i 

Let C supj.>o Gl_^/r'^Zi- Now implies that Gu/Gk-i < 1 + AG. From 
this and (|5.6p we get 

G\_i Gk-i ^ ^ (^Gfc_i Gk 

Combining all of the above and (|5.5p we obtain 

(1 + AG) 



< 



- Ae 



□ 



We now prove Theorem [T] 



Proof Let /3 = supo<f<T &(||e„(t)||oo) (by (A2), this is finite as ||e„(t)||oo is contin- 
uous along an orbit). Then (A5) becomes for u{t), t e [0,T], and all a £ , 

f < Pd- 

We now write the integral form of this and get (recall that |n*(Q;)| < ^fN): 
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F\R,T) =11 f{a,t)n{a)dadt 



T 

2/ 



< T\dC*{r)\-sjN I / r{a,t)dadt 

JdC"(R) 



< T\dC* {R)\Vn ■ (3 [ [ d{a,t)dadt 

Jo JdC'iR) 



We introduce the constant 

(5.8) LON = ^3/22^, 

and then \dC*{R)\y/N — a;Ari?^^^As dC*{r) are disjoint, we now have 

1 '^F^{r,T) 



(5.9) D{R,T)> 



r—l 



Combining that and (|5.4p . we see that for all R E Jt, 

1 ^ F^ir,T) 



F{R,T) > 



r— 1 



^ reJT,r<R-l. 

As u is not stationary, there exists Ri large enough so that D{Ri,T) > e for 
some £ > 0. If Jt is finite, the claim is in both cases proved. If not, by (|5.4I) there is 
Ro G Jt, Ro > Ri and then F{Ro,T) > e > 0. Now if Rq = ro < ri < r2 < ... are 
all elements of Jt greater or equal than Rq, we can define a recurrent sequence 

(5.11) Go = e>0, 

(5.12) Gk = , , arj. ^ = Gk-i + „ j^I^ ■ 

From this and (|5.10p we easily deduce by induction that for all fc > 0, < 
i^(rfe,T). However, by (A2), t ^ ||/u(t)||oo is continuous, thus bounded on [0,r] 
by some constant G. By definition then we have for all r > 0, thus F{r, T) < 
Cuj^r^^^ . We deduce that Gk/r^~^ is bounded. The claim now follows from 
Lemma S) □ 

If the dimension iV = 1 (as already shown in [21]), the theorem[l]can be rephrased 
as follows: given any non-stationary u E X and any T > 0, there exists Rq such 
that for aU R > Rq, 

(5.13) E{R,T) < E{R,0). 

Note that the choice of Rq is typically non-uniform: it depends on T > and 
the choice oi u £ X. 

In dimension N = 2, the set of all R for which (|5.13p holds has density, as the 
condition that X^r-e/r ^^''^ finite is equivalent to 

K 



\n{u,T) n{i,. .,K}\ = o 



logK 
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Remark 1. In dimensions N > 3, the rationale above does not bring any new 
insight, as X^rGN ^/''^^^ '■^ always finite. 

Corollary 1. // (^,1^3) is a Lattice EDS in dimensions N = 1 or 2, it has no 
periodic orbits of period T > 0. 

Proof For periodic orbits of period T > 0, E{R, T) = E{R, 0) must hold for all R, 
which is in contradiction with Theorem [TJ □ 

6. The flux and dissipation bounds 

In previous section, we investigated constraints on the energy of a cube C*{R) 
at a time T > 0. Here we will use similar rationale to develop bounds on the total 
flux F(i?, T) through a boundary dC*{R) of the cube up to a time T > (see ([O)) 
for the definition of F[R^ T)). In this Section we will assume that a Lattice EDS is 
bounded, with the constant /? such that (|3.7p holds. We fix an orbit u{t), and use 
the notation eo = ||e„(o)||oo. 

We first state the flux bounds which will be proved in the following. Let T > 0, 
and N as usual the dimension of Lattice EDS. As we are typically interested in the 
asymptotic behavior of our Lattice EDS, we develop the bounds valid for large T 
as compared to R. 

(i) : UN =1, then 

(6.1) lim sup -^F{R,T) < 2y^. 

(ii) : If A^ = 2, then 

(6.2) lim sup ^-^^F{R, T) < 96\/2 • /3. 

T->oo J 

(iii) : If iV > 3, then 



(6.3) lim sup -F{R,T) = n'^/^2^n-2 . 



T-i-oo T 



(3R 



N-2 



Let us first discuss the significance of the bounds above. A simple a-priori anal- 
ysis can yield only the bound F{R,T) = 0{TR^^^), i.e. that for bounded Lattice 
EDS, the maximal flux through a boundary of a cube C]j is at most proportional 
to time and the size of the cube. The bounds above are much better. In low di- 
mensions N = 1,2, we see that the flux can not be proportional to the total time 
spent. We will see later that because of that, our Lattice EDS will be for almost all 
times arbitrarily close to equilibria. In the dimension > 3, we see that the flux 
can grow proportionally with time, but only proportionally with R^~^ . This will 
mean that our dynamics may have persistent non-equilibrium behavior only on the 
sets of the size ^ R^~^, i.e. two dimensions less than the space dimension. 

We will prove the flux bounds by showing that F{r,T) < Gr, where Gr is a 
recurrent sequence defined for integer r, satisfying 

(6.4) Gr+i - Gr-Xr^'-'+e-^, 

(6.5) G,./r^-i is bounded. 

As we will see later, here A ^ eo, £ ~ l/(/3T). The analysis of the recurrent 
sequence for different A^ is similar to the proof of Lemma 21 but significantly more 
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technical. Therefore the proof of the following two Lemmas is given in the Appen- 
dix. 

Lemma 5. Assume Gr, r > 1 satisfies (6^, 116. 5\) for some X,e > 0. 
(i): IfN^l, then 

(6.6) Gr < 



(ii): If N>2, th. 



en 



(6.7) < (7V-1) ( 1 + - ) — + \j-r 



N-1 



As we will see in the Appendix, the bound (|6.6I) in the case N = lis sharp, and 

(|6.7|) in the cases iV > 3 is sharp up to a constant. We will be interested in the 

asymptotic case e — (i.e. T — )■ oo), so in the case N = 2 we can say something 
more. 

Lemma 6. Assume G^, r > 1 satisfies i6.4\ l, i6.5]) . Assume also N — 2 and 
ry/Xe < 1/2. Then 

(6.8) Gr < ^'^ 



-elog{2r^Xe) ' 

We will use it to prove the following general flux bounds: 
(i): If iV = 1, then 



(6.9) F{R,T)<2^p^. 
(ii): If iV = 2 and MR^eo < W2/3T, then 

(6.10) F{R,T)< ^'^ 



log(w2/3r/(64eoi?2))- 
(iii): If TV > 2, then 

2 ^ 



(6.11) F{R,T) < (N-l) i^l + -j c^Ari?^-2^T + 2^V^A, • N-R^'-'y/^. 

Theorem 2. Assume a bounded Lattice EDS is given, with the constants P and eg 
as above. Then for allT > 0, R>1, fOZ!)] and fgjT]] hold. 



Proof. The first step is writing again the energy balance equation. As |C*(i?)| 

.(0) 



2^i?^ and eo = ||e„(o)|U, we have E{R,0) < 2^ R^ e^. Inserting that and e > 



in (15. 4p we get 
(6.12) F{R,T)>~2^R^eo + D{R,T). 

As in the proof of Theorem [1] (j5.9p holds. Inserting it into (|6.12p . we get 

F(i?,r) > -2^R^eo + ^J2^'^''^^ 



unPT ^ r 

r—l 



N-l 



R-1 



> -2'^eoNj2r^-' + ^Y: 



r=l 



DISSIPATIVE SYSTEMS ON LATTICES 



17 



where we used 

Jo 



Now we can set 

s = 1/{UNI3T)- 
We define the recurrent sequence for some given Gi € 

a 



(6.13) Gn = -Xj2r''-'+sJ2 , 

r—l r—1 

From this, it is easy to see that for all r > 1 we have 

Now, if for some tq, Grg = F{ro,T), inductively we have that for all r > ro, 
Gr < F{r,T). We can also see that F{r,T)/r^-^ is bounded, as the flux ||/||oo is 
uniformly bounded on [0,r] along the semiorbit u{t) by (A2). Thus we have 

F(r,T)<wjvr'^"i sup ||/„(t)||oo. 

tG[0,T] 

We conclude that sup^>o Gr/r^^^ < oo. We flnally apply Lemmas [S] and [SI and 
insert the definitions of A,£ into (|6.6p . (|6.7p and (|6.8p respectively. □ 

The asymptotic behavior for large T as compared to R now straightforward: 

Corollary 2. Assume a hounded Lattice EDS is given, with the constants j3 and 
eg as above. Then \6. h6.^) and \6.S\l hold. 

Proof. We insert into ((O)) . (pl^ and (pHt the definition ([5?S)) of wat, and in the 
case TV > 3 also (1 + < 2^-2_ □ 

We can now deduce the following dissipation bounds. 

Corollary 3. Assume a bounded Lattice EDS is given, with the constants /? and 
eg as above. Then for any integer R > 1 
(i): IfN = l, then 

i-T 



(6.14) lim sup I I d{a,t)dadt — 0{\//3eo). 

T-i-oo vT Jo Jc*(R) 

(ii) : If N = 2, then 

(6.15) lim sup / / d{a,t)dadt = 0{^). 

T^oo J Jo Jc'{R) 

(iii) : IfN>3, then 

(6.16) lim sup / / d{a,t)dadt^O{PR^-^). 

T->oo T Jo Jc*{R) 

Remark 2. The notation F{x) < 0{G{x)) means that there exists a constant C > 
such that for allx, F{x) < CG{x) (and not, as is more common, |i^(a;)| < CG{x)). 
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Proof. Let = 1. The claim can be written as 

lim sup -^D{R,T) = 0(7^) 
T->-oo yT 




(|6.3p instead. 



□ 



These bounds on dissipation have important interpretations. In dimensions 1, 2, 
we see that for an average < > 0, the total dissipation on any chosen cube C*{R) 
at the time t is arbitrarily close to 0. We can write it shortly 



We will see in the next section that it means that any semiorbit of a bounded 
Lattice EDS stays for almost all times arbitrarily close to the set of equilibria. 

In dimensions iV > 3, this conclusion does not necessarily hold. We see, how- 
ever, that the dynamics which can persist on the set of times of positive density is 
concentrated on a set of two dimensions smaller than the dimension of the physi- 
cal lattice space. For example, in dimension N = 3, the non-stationary dynamics 
of dissipative lattice dynamics which can persist is at most one-dimensional, thus 
essentially concentrated on a finite set of lattice points. 



Now we deduce dynamical consequences of the bounds on flux and dissipation 
deduced in the previous section. We focus for now to dimensions N = 1,2. We will 
assume that our Lattice EDS is bounded. 

Also we assume that the Lattice EDS semiflow {X, ip) is given on a compact 
space X. This is not overly restricting: we have shown that, for example in the 
case of the FK model, all orbits of bounded width are relatively compact in the 
quotient space X = Y j where with ~ we identify all configurations which are 
integer translates. The topology here is the induced product topology, i.e. the 
topology of pointwise convergence on lattice points. 

We can thus study the continuos semiflows 93 on compact metrizable spaces X, 
with the structure of a bounded Lattice EDS. In this setting, we are interested on 
the structure of w-limit sets of any m G X, i.e. the set of limit points of u(t) as 
t — > 00 in the chosen topology on X. 

We first show that the average time a flow spends outside of any neighbourhood 
of stationary points is 0. Denote by E CX the set of all equilibria (or stationary 
points) of the semifiow 1^, that means all points u ^ X such that for all t > 0, 
u{t) = (p{u, t) = u.liU C X is any set, then 1^^ : X M is as usual its characteristic 
function, lu{v) = 1 if u € otherwise lu{v) — 0. 



(6.17) 



lim sup -D{R,T) = 0. 

T^oo T 



7. Convergence to equilibria in dimensions 1,2 
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Proposition 1. Assume {X, ip) is a bounded Lattice EDS, assume X is compact, 
and let N = 1 or 2. Let U be any open neighbourhood of the set of stationary points 
£ . Then 

(7.1) lim 1 r lu{u{t))dt^l. 

T^oo T Jq 

In particular, S is non-empty. 

Proof. As X is compact, tlie complement Wis also compact. Choose u G W. As 
u is not stationary, by (A3) there exists R{u) large enough and e{u) > such that 

/ du{a)da > 2e{u) > 0. 

Jc-{Riu)) 

By (A2), that means continuity of u o?„(q;), there exists an open neighborhood 
V{u) oi u E such that for each v E V(u), 

dv{a)da > e{u). 

The family V(u), u eW^ is an open cover of the compact set W^, thus it has a finite 
subcover V(wi), V{uk). Let s — min{e(ui), e(?ife)}, R ~ max{i?(ui), R{uk)}. 
Then for any u G W^, 

/ du{a)da > e > 0. 
Jc"{R) 

We now have 

1 1 f^i 



lim sup — / luc{u{t))dt < lim sup — / -l^c{u{t)) / du{t){a)dadt 

T^oc J- Jo T->oo J Jo £ JC*{R) 

< -\im sup ^D{R,T), 

which is by (|6.17p equal to 0. □ 

The LaSalle principle for gradient systems states that the w-limit set of each 
point consists of stationary points. Our weaker version of this for Lattice EDS is 
as follows. 

Proposition 2. // {X, ip) is a bounded Lattice EDS, assume u[t) for some u E X 
is relatively compact, and let N — 1 or 2. Then u){u) contains a stationary point. 

Proof. We assume without loss of generality that X is equal to the compact closure 
of u{t), thus compact. By Proposition[l] u{t) visits any neighborhood off infinitely 
many times. By compactness, u{t) than has a limit point in £. □ 

We recall the notion of uniform recurrence, as the strongest notion of recurrence 
which occurs in general dynamical systems. A point x € X is uniformly recurrent 
with respect to a continuous semiflow </? on a metrizable space X, if for each open 
neighbourhood U of x, the set of return times r = {t > | f{t,x) £ C/} satisfies 
the following: 

R: There exists n > 1 andti,...,i„ so that [0,+oo) = UjJ^^^jt > 0,t + tk e U}. 
It is easy to show that Proposition [1] implies that the set of uniformly recurrent 
points of a bounded Lattice EDS coincides with the set of equilibria £ (see e.g. 
Lemma 5.6 in [IS])- We also recall the relationship between minimal sets and 
uniformly recurrent points (see Furstenberg |12] . Lemma 1.14 and Theorem 1.15): 
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Theorem 3. Given a continuous semiflow on a compact metrizahle space, every 
minimal closed positively invariant set consists of uniformly recurrent points. 

We can now summarize the key dynamical systems properties of bounded Lattice 
EDS on compact sets in dimensions 1,2. 

Theorem 4. Assume {X, (p) is a 1 or 2-dimensional bounded Lattice EDS, and 
assume X is compact. Then 

(i) : {X,(p) has no periodic orbits of period T > 0; 

(ii) : The uj-limit set of each u X contains an equilibrium; 

(iii) : The only uniformly recurrent points are equilibria; 

(iv) : The only minimal, closed positively invariant sets consist of single equi- 
libria; 

(v) : All invariant (Borel probability) measures are supported on the set of 
equilibria. 

Proof. We proved (i) in Corollary [TJ and (ii) in Proposition [5J The definition 
of uniform recurrence and (|7.1I) imply (iii); (iv) follows from (iii) by Theorem [S] 
Relation (|7.ip and the Birkhoff ergodic theorem imply (v). □ 

Discussion of dimensions N > 3. A natural question to ask is whether some 
claims of Theorem [3] can be extended to dimensions TV > 3. We have already noted 
in Remark [T] that the constraints to the flux transfer in dimensions iV > 3 are 
qualitatively different, as the sum X^rGN ^/^^^^ is finite. 

Indeed, in dimensions > 3 one can construct energy increasing orbits and 
periodic orbits, thus contradicting all the claims of Theorem |4l The examples in 
the continuous space case have been constructed in |13l , and can be adapted to the 
discrete space case. 

We believe, however, that the asymptotic bound on dissipation D{R,T) ~ 
Ti?^-2 as T » i? in dimensions A > 3 can imply interesting conclusions for par- 
ticular equations. We suggest in Section [9] an application to coarsening dynamics 
in bistable potential. 

Structure of the w-limit set in dimensions 1, 2. Another important question 
is whether one can describe in some detail the structure of w-limits set of compact, 
bounded Lattice EDS in low dimensions. One can show that even in dimension 
1, for compact, bounded Lattice EDS, the w-limit set can contain non-stationary 
points, so Lattice EDS indeed differ from the gradient systems and lattices on 
finite domains. As in the continuous case (|14j). our example is inspired by the 
coarsening dynamics of the real Ginzburg-Landau equation studied by Eckmann 
and Rougemont, and is outlined in section |9l 

8. Bounds on relaxation times 

We can use the flux bounds ()6.9p . ()6.10p and ()6.1ip to calculate upper bounds 
for an arbitrary orbit to relax e-close to the set of equilibria f on a cube of size r. 
Given a bounded Lattice EDS, let U^^r be the set of all u e AT such that for all 
a e C*(r), d(a) < e. Then by (A3), 

e>0, r>l 

Let te^rbe the first time such that u{t) e U^^r- Let /3, eg = ||e.u(o)||oo be as in the 
Section [9l We will show that: 
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(i) : For iV = 1, i,,, < eoP/ie^r'^), 

(ii) : For = 2, U.^ < exp(c2/3/(£r2)), 

(iii) : For iV > 3, 

eo/3 



te.r < 



whenever the denominator is > 0. 

Here cjv are constants depending only on the dimension A^, and are evaluated 
below. 

If u(t) ^ Ue^r on [0,T], the energy balance equality, e > and \C*{r)\ = 2^r^ 
imply, after integrating over [0,T]: 

(8.1) F{r, T) + eo2^r^ > D{r, T) > e2^r^T. 

The case = 1. It is easy to check that, if AT - B^/T - C > for some 
constants A,B,C > 0, then 

(8.2) Vt < B/A+ y^CjA. 
We combine that with (|6.9p and (|8.ip and get 



The case N — 2. Again we first consider a general inequality AT/(logT— _B) 
C > DT, for some constants A, B,C,D > 0, which implies 



thus 



, A B\ C 
T < max < exp — — - + — , — =- 



2D 2 J ' 2D 

From (PUT)) and (gl]) in the case B = log(64eorV(^^2/3)) > -4=^ Meor"^ > uj2l3, 
we have 

(32eor^ ( 2u:2^\ eo 

te^r < max <^ — exp — , — 

I UJ2P \ er-' ) 2£ 

(the case B < Q results with only a small order correction, and can be evaluated 
easily). We see that in dimension 2, we have two relaxation timescales: the fast re- 
laxation for e of similar magnitude as the locally available energy eo, and potentially 
exponentially long relaxation time for e ^ cq. 

The case > 3. Similarly as above, we use (|6.1ip . (|8.1I) and (|8.2I) and get 

NujNeoP ( eoy/P \ 



r2 (Ara;Ar/3/(4r2) - ef \r {NujNl3/{4r^^ - -^^'^ 



e] 



if e > NujN(3/{4r'^). If e < N lo n (3 / [Ar"^) , the orbit does not necessarily ever enter 
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9. Application example: Coarsening in bistable potential 

9.1. Introduction and setting. Eckmann and Rougemont studied in [7] the real 
Ginzburg-Landau equation 

(9.1) Ut = Uxx + u-u^ 

(a continuous analogue of an elastic Id lattice in a polynomial potential), also as an 
extended system, i.e. considering time evolution of configurations u : R — > [—1,1] 
not necessarily vanishing at infinity. They were interested in dynamics of multikink 
solutions: evolution of initial configurations u(0) which "jump" between minima ±1 
infinitely many times. They then have shown that for selected initial conditions, 
for each a; G M, u{x) jumps infinitely many times between ±1. They also argued 
that the resulting dynamics of lengths of intervals for which a configuration u lies 
in the same minimum approximates well the abstract Bray and Derrida coarsening 
dynamics model [6]. Rougemont in |19j studied behavior and disappearance of 
selected initial "droplet" conditions for the analogue of (|9.ip in dimension 2. 

Here we develop a lattice equivalent of this model and behavior, in somewhat 
more general setting. In particular, we will also discuss second-order damped dy- 
namics in dimensions 1,2. 

Assume L : R x M — > M is a smooth function (at least C^) satisfying (LI), (L2) 
from Section m and in addition the following (satisfied e.g. by the standard Id FK 
model) : 

(L3): The pair (0, 0) is the strict global minimum of L; 
(L4): The partial derivative L12 > 0; 

(L5): L(x,y) = L{-~x,-y). 

We will consider the dynamics on the A'^-dimensional, scalar valued lattice, with 
the formal energy associated to u : 7L^ — >■ M: 

N 

We thus discuss the dynamics 

N 

(9.2) A • dttu{a) + dtu{a) = - ^(L2(u(a - ej),u{a) + Li{u{a),u{a + Sj)), 

where A > 0. We will distinguish below the case A = (gradient) and A > 
(damped). Both cases belong to a class of generalized FK models, discussed in 
Section [H with the standard Id FK model being a special case. In the damped 
case, we will further require that A is not large enough, so that the partial order 
introduced by Baesens and MacKay [3], [1] holds. Let 

B= max {L22{x,y) + Lii{y,z)). 

0<x,y,z<l 

We will assume the following (trivially holding for A = 0) overdamped condition: 
(L6): Assume 4XB < 1. 
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9.2. Dynamics of ()9.4p as a bounded Lattice EDS.. Wc will be interested in 
solutions < u < 1. Because of periodicity of L and (L3), Wo,wi E X\ as defined 
below are stable equilibria: 

A = 0: wq = 0,wi=0, 

A > 0: TWO = (0,0),u;i = (1,0), 

where in the case A > we consider as usually at a single lattice point the evolu- 
tion of w{a) — {u{a) , dtu{a)) . Analogously to Eckmann and Rougemont, we will 
be interested in solutions of (j9.2p "jumping" at each lattice point infinitely times 
between wq and Wi . Thus the following choice of the state space is appropriate: in 
the case A = 0, we will assume that for all a E , 

(9.3) 0<u(a)<l, 

and in the case A > 0, we will in addition require that the derivative v{a) = dtu{a) 
is not too large: 

(9.4) < 2A- u(a) +M(a) < 1. 
The state spaces will thus be as follows: 

• A = 0: Xo = {u:Z'^^R satisfies ([931) }. 

• A > 0: = {{u, u) : X ^ M satisfies (lO)) . (jOll } . 

Clearly by Tychonoff theorem, for all A > 0, X\ is a compact set in the induced 
product topology (i.e. topology of convergence at each lattice point). By combining 
results from Section |4] and the order-preserving property of the solution of (|9.2|) , 
we see that the theory developed in this paper applies to (|9.2[) . 

Proposition 3. Assume A > and L is satisfying (L1)-(L6). Then the solution 
of \9.2\l generates a bounded Lattice EDS on Xx, which is compact in the induced 
product topology. 

Proof. We first show that X\ is invariant. Note that the condition (L4) implies 
that the lattice with the dynamics (|9.2I) is cooperative (see relation (2) in [T). In 
the case A = 0, the ordering on M^" is defined with u < u*, if for all a E Z^, 

u{a) < u*{a). 

By [3], Proposition 2.1, the dynamics (|9.2[) preserves ordering, thus wq < u{Q) < wi 
implies that for all t > 0, wq < u{t) < wi, so Xq is invariant. 

Let A > 0. The partial order on (M x M) introduced in [3] is now (u, v) < 
(m*,?;*), if for all a e Z^, 

u{a) < u*{a), 
2Xv{a) + u{a) < 2Xv*{a) + u*{a). 

We set as usual v{a) = dtu{a). By 4 , Proposition 4.2 and (L6), the equation 
(|9.2p preserves the partial order above. The conditions (|9.3|) . (|9.4p can be read as 
Wq < {u{0) , dtu{a)) < wi, thus by the order preserving property, for all t > 0, 
Wq < {u{t) , dtu{t)) < wi. We conclude that Xx is invariant. 

We have shown in Section 2] that (19. 2[) generates a Lattice EDS. By compactness 
and (|9.3p . the energy 

N 

eu(a) = L{u{a), u{a + Sj)) 
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is bounded, thus we have a bounded Lattice EDS. □ 

9.3. The set of equilibria. The set £ of equihbria of (|9.2I) corresponds in the case 
A = to the solutions of 

N 

^(L2(u(a — Ej), u{a) + Li{u{a),u[a + Sj)) = 0, 

and in the case A > to the same set with v ~ Stu = 0. For the moment, we 
consider the case = 1, and assume (L4) is somewhat stronger, i.e. there is 
(5 > so that Li2 > S. By the well-known correspondence between the equilibria of 
generalized Id FK models and orbits of area-preserving, positive twist maps, there 
exists a map on the annulus / : xR defined with (u,p) {u* ,p*), p* ~ L2{u, u*), 
p = -Li{u, u*) (see e.g. [T7], Section 9.3.b). 

We will be interested in functions L which generate twist maps with a phase 
portrait as the classical Id pendulum (see Example [T] below) , namely the cases 
where stable and unstable manifolds of the equilibria (0,0), (1,0) of the twist map 
consist of their heteroclinic connections. 

Example 1. We choose a smooth potential : M — > M, periodic with period 1, 
symmetric (V {x) = V{—x) ), with unique maximum per period at x = 0; for example 
V{x) = sin(27ra;). Then we consider the continuous flow ip of {u, Stu) of the second 
order equation 

duu{t) = V'{u{t)) 

on the annulus x R, and /et /v : S"'^ x R — S"*^ x R he its time-one map. As fv 
is an area-preserving twist map, one can define its generating function Ly which 
satisfies (L1)-(L6) (see also [17 , Section 9.3.h). 

This motivates us to define the set of stable or asymptotically stable equilibria 
S <Z £ (ZX\ containing wq^wi^ and all v G for which there exists a vector 
e £ {±£1, zbe^v} so that 

(9.5) lim T'^v e {wo,wi} 

n^oo 

(in the product topology). In the Example [U the set S is connected, so, as we will 
see, a solution of (|9.2I) can "slide" between the equilibria wq and wi. 

9.4. Translationally invariant ergodic measures. Rather than discussing as- 
ymptotic behavior for a particular initial condition in X\ , we will prove existence of 
coarsening dynamics for almost every initial condition with respect to a particular 
Borel probability measure on X\ . We now specify properties of such measures and 
construct some examples. We will always consider Borel probability measures on 
X\ with respect to the induced product topology on X\. 

The group naturally acts on Xx as the group of translations T^^ for /3 e Z^: 

For A = 0, (r^u)(a) = it(a-f /?), 

For A > Q, {Tf^{u,v)){a) = {u,v){a + P). 

We say that a Borel probability measure /i on X\ is translationally invariant, 
if for each /3 € and each measurable A C Xx, fi{T^{A)) = p,{A) (clearly it is 
sufficient to show this only for generators T^^ , T^"). A set A C Xx is trans- 
lationally invariant, if for each (3 e Z^, T'^(A) — A. A translationally invariant 
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Borel probability measure /i is ergodic, if all measurable translationally invariant 
sets have measure or 1. 

We can also naturally define the reflection operator R on X\ , which " swaps" the 
stable equilibria and 1: 

For A = 0, {Ru){a) = 1 - u{a), 

For A > 0, {R{u,v)){a) = {l-u{a),-v{a)). 

Note that (LI) and (L5) imply that the equation ()9.2|) is i?-invariant. The first two 
required properties of Borel probability measures fi on X\ are the following: 

(Ml): /i is translationally invariant and ergodic, 
(M2): /i is i?-invariant. 

There are many Borel probability measures on X\ satisfying (M1),(M2); we 
construct only one, canonical example fi^^, A > 0. 

Let A = 0. Let f be the BernouUi measure on {0, 1}, J^({0}) = i^Cjl}) = 1/2, 
and let = i'^ (defined in the standard way as a product of countably many 
probability measures). This can be understood as randomly and independently 
putting each lattice point in one of the minima 0, 1 with the same probability 1/2. 
If A > 0, we do the same, and set the initial velocity to be almost surely 0. More 
precisely, ^j-x = HqxS^q^, where (5{o} is the atomic measure such that (5{o}({wo}) — 1, 
Mo = 0. 

Now fi^ are by definition translationally invariant and i?-invariant. It is easy to 
show that fix is ergodic (e.g. by applying i5^. Theorem 1.3). 

9.5. Existence of coarsening. In addition to (Ml), (M2), we also assume the 
following 

(M3): For fi-a.e. u e Xx, w(m) n (f - 5) = 0. 

Here U!{u) is the w-limit set with respect to dynamics (|9.2p . In dimension = 1, 
this means that L generates a twist map with a phase portrait as in Example [TJ 
and that fi-a.e. u does not lie on a stable manifold (with respect to the dynamics of 
(|9.2p ') of some of spatially periodic or quasiperiodic equilibria different from iuq, u^i. 
As discussed also in [7, in the PDF case, this is a reasonable assumption but a 
technically difficult claim. In [22] we will show the following: 

Claim 1. The measures fix constructed above satisfy (MS) in dimensions N = 1,2 
for the family of functions Ly constructed in Example [7J 

Here we bypass this difficulty, and prove the following existence of coarsening 
dynamics in the general, abstract case, by assuming (M3): 

Theorem 5. Assume X > 0, L is a function satisfying (L1)-(L6) , and fi is a 
Borel probability measure on Xx satisfying (M1)-(M3). Then for fi-a.e. u € Xx, 
its uj-limit set with respect to i9.S\) contains both stable equilibria wo,Wi. 

Proof. Let Aq be the set of all u £ Xx so that uj{u) contains either wq or an 
equilibrium v asymptotic to wq in the sense of (19. 5p . and let Ai C Xx be an 
analogous set with respect to wi. By Proposition|3l the dynamics of (19. 2p generates 
a bounded Lattice EDS semiflow on a compact set Xx, thus Proposition [7] applies. 
So for each u G Xx, its w-limit set contains an equilibrium. Now because of (MB), 



li{AoUAi) = 1. 
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However i?-invariance of ()9.2|) and (M2) easily imply that ^{Ao) = ^{Ai). As 
Ao,Ai are translationally invariant, by ergodicity, fi{Ao) , fJ.{Ai) e {0, 1}. We con- 
clude that /i(Ao) = fJ-iAi) = 1. Finally, by the translational invariance of Aq^Ai, 
relation (|9.5p . and closedness of w- limit sets, we deduce that each u ^ Aq must con- 
tain wq in its w-limit set, respectively each u £ Ai must contain wi in its w-limit 
set. Now ^0 n Ai is the required subset of X\ of full measure. □ 

The interpretation of Theorem [5] is that for /i-a.e. initial condition, the size of 
"droplets" of lattice points in the same equilibrium or 1 grows with time, but a 
single lattice point changes its position infinitely many times. By Corollary [3l the 
average time a lattice point spends "moving" between wq and wi is 0. One can 
also use the bounds on relaxation times in Section [S] to calculate upper bounds on 
the maximal "lifetime" of droplets of certain size. 

9.6. Non-equilibrium point in a w-limit set. As was mentioned in the intro- 
duction, the dynamics of Lattice EDS topologically differs from gradient systems, 
as we can construct initial conditions u so that uj{u) contains a non-equilibrium 
point. This is possible even for bounded Lattice EDS on compact sets in dimension 
= 1. We consider the equation (|9.2|) . let A = 0, dimension A^ = 1, and assume L 
is of the type Ly constructed in ExamplelU so Ly satisfies (L1)-(L6). As discussed 
in more detail in [14], Example 5.7, in the continuous-space case, we consider an 
initial condition 

.X / b2n < |a| < hn+l 

u(a) = < , , I , 

^ ^ \ 1 b2n+l < \a\ < 62n+2, 

where = 6o < < &2 < •■• is an increasing sequence of integers, satisfying also 
bn+i ^ bn- As explained in |14| . the annihilation of kinks in the coarsening sequence 
always happens at the origin a = 0, so for a sequence 6„ increasing rapidly enough, 
w-limit set of u in Xq contains both wo,wi but also non-equilibrium, heteroclinic 
orbits of (|9.2p connecting them. A detailed discussion and a proof will appear in 
[22]. 



10. Appendix: Proofs of Lemmas 1316] 
Here we study the recurrent sequence 

(10.1) Gr+i - Gr - Ar^-i 



such that Gr/r^ ^ is bounded. We substitute gr = Gr/r^ ^, and can study it as 
orbits of a discrete dynamical system G{r,g) — (r*,(7*), 

(10.2) r* = r + 1. 



(10.3) g* = ig-X + eg') 



, 7V-1 
2n ' " " 



1 



where A,e > are parameters. G is defined on [l,-|-oo) x R — )-[l,-|-C!o) x R. We 
want to give sufficient and sharp conditions on (r, g) so that the second variable is 
bounded along the orbit of G. We first note that for r = oo, G has two fixed points 
±y/X/e. We will see that the fixed point +^/X/e is asymptotically a saddle point, 
and that the necessary and sufficient condition for 5 > to be bounded along the 
orbit of G is that (r,g) is below the stable manifold of +^A/e. More precisely: 

Lemma 7. There exist a continuos curve gs : [l,+oo) x MU {+00} such that: 
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(i) : Ifg = gs{r), then lim„^oo G"{r,g) = (+00, y^Xjs). 

(ii) : IfO<g< gs{r), then lini„^oo G'"(r, .g) = {+00, ~y^^). 

(iii) : Ifg > gs{r), then lini„^oo G'^{r,g) = (+cx),+oo). 

(iv) : The curve gs{r) is non-increasing. If N > 2, gs{r) is strictly decreasing. 

Note that we did not exclude the possibihty that the unstable manifold gs di- 
verges for some ro > 1 (and is then formally = +00 on [l,ro]. We will see later 
that this is not possible. 

Proof. In order to apply standard results on invariant manifolds of fixed points, we 
need to choose the right substitution and compactify (I10.2[) . (I10.3P at r = +00. We 
set the bijective map [0, 1] — > [1, +00], s i-> r with 

1 



1 



If i{s,g) 
R^[0,1] X R such that G 
is given with (s, g) 1 — > (s* 



{r,g), i ^{r,g) — {s,g), we now analyze the map H : [0,1] x 



^ o H o i. Then it is easy to check that the map H 



2-s' 



9 



(.9 - A + eg2) 



1 



N-l 



Clearly H is on [0, 1] x R and that it can be extended to a map on a 
neighbourhood of [0, 1] x R. As s* > s on [0, 1), the map H has two fixed points 
{s,g) — (1, ii/A/e). We calculate the differentials at fixed points, and get 

1 
1 



DH{l,y/X/e) = 



Clearly the fixed point (1, — i/A/e) is a sink: the s = 1 direction is attracting, as 
the direction s — I corresponds to the eigenvalue (0, 1) and the eigenvector 1 — 2\/Ae 
and is thus attracting, and the other direction is attracting as s i-> s* is increasing. 
We consider now the fixed point (1, y^X/e) with eigenvalues 1,1 + 2y/Xe. We apply 
the central manifold theorem, and construct an unstable and a central manifold. 
Clearly the unstable manifold is s = 1. Assume a local central manifold is given by 
{s,gs{s)), : [1 — (5, 1] — >■ R (because of the expression for s s*, we can without 
loss of generality assume that the central manifold is locally parametrized by s). 
One can check that the central manifold is smooth (as H is smooth in the neighbor- 
hood of (1, ^J\/e)), locally unique, and that it can be extended to a global smooth 
central manifold (s, (7s(s)), gs '■ [soi 1] — > R, < sq < 1 by extending it with values 
of H^^ for which 5 > 0. Also if sq > 0, then lim^^^+ gs{s) — +00. Furthermore, as 

s s* is strictly increasing, for each {s,gs{s)), lim„^oo ^^"(s, (7s(s)) = (1, yJX/e), 
so the central manifold is asymptotically stable. 

We now show that gs{s) is decreasing. From the direction of the central eigen- 
vector at the fixed point (1, y^X/e), we deduce that the function gs{s) is eventually 
decreasing (strictly decreasing if A'' > 2). However, one can easily check that 
if for any (si,5i), (32,32) = H{si,gi), (53,53) = i?(s2,52), 92 > gi > implies 
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53 > 32 -(and 53 > 52 if s < 1 and TV > 2). From this we see that gs{s) is decreasing, 
and in the case N > 2 strictly decreasing. 

Consider now the negative image of the unstable manifold, i.e. the maximal set 
(s,5'~(s)) satisfying H{s,g~{s)) = {s*,gs{s)) and gj{s) < 0. Then we can divide 
the state space [0, 1] x M of the function H into five sets: the sets of {s,g) such that 
A = {g> ssis)}, Ws = {g = gs{s)}, B = >g> gjis)}, = {g ^ gj{s)} 

and C = {g^ (s) > g}. One can now summarize the discussion above, and show the 
following: H{A) C A, H{Ws) C W,, H{B) C B, H{W~) C Ws, H{C) C A, and 
(1, y/X/e) G Ws, (1, —y/X/e) S B. As the dynamics of H asymptotically coincides 
with the dynamics along the invariant set s ~ 1, we get the following: 

is,g) G A(J C ^ Mm H"{s,g) ^ {1, +00), 

n^oo 

is,g) e WsU W- ^ lim H'\s,g) = (1, ^X/^), 

n— )-oo 

e lim i/"(s,,g) = (l,-VA7i). 

n— )-oc 

By substituting back r instead of s and writing gs{r) — gs{s{r)), these conclu- 
sions imply all claims of the Lemma. □ 

We now prove Lemma [5j We now see that it is necessary and sufficient to 
give upper bounds on Gr on the stable manifold, that means of the form Gr = 
gs{r)r^~^ ■ The reason is that by Lemma[71 if G,. > gs{r)r^~^ , limr^ooGr/r^^^ = 
+00. Also if < < gsir)r^-\ lim^^oo GJr^-^ ^ ±y/X/£- 

In the case iV = 1, the stable manifold can be expressed exactly and is gs{r) — 
y/ A/e, which implies (|6.6p . 

We deduce the case N > 2 from the fact that gs{r) is strictly decreasing, thus 
gs{r + I) < gsir). That means 



gs{r) - X + egsirf < gsir) 
1 + y/x for X > we get 



If Cr — (i-ti)^ ^ — 1, by solving the quadratic inequality and applying y/1 + x < 



Gr+Cr^l+4SX/G^ Cr /A 

(10.4) gsir) < < - + Vl- 



However 

(10.5) Gr < ' ' 



N-2 



Combining (ITU^ . (ITU3)) and Gr < 5s(r)r^"i (^TT)) . Also gsir) is finite for 

all r e [l,+oo). One can also show by constructing appropriate examples, that 
(|6.7p is up to a constant the best upper bound in cases iV > 3. 

The case N = 2 and Lemma [6] is, however, more subtle. Fix N = 2, and assume 
Gr = gsir)r. As Gr > ry/xje > 0, we divide (|10.ip by G^ and get 
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Inserting Gr > X/e in (llO.ip . we get Gr+i > Gr- However, as gs{r) is 
decreasing, Gr+i/Gr = gs{r + l)/gs{r) ■ (1 + l/r) < (1 + 1/r) < 2. We now have 

/ _^ Gr4-^ 1 Gi — LI Gr Gr-i-^ Gr /I 1 

(10-7) ^ - 7^ = < 2 =27^- 



G^ GrGr-\-1 \Gr Gr-\-i_ 

We now fix an integer k > 2, and sum ()10.6p for s = r, fcr — 1. Applying (jl0.7p 
and Gr+s > Gr we fiave 

kr—l kr- 



11 A X - e X - 1 



1 

Gj- Gr Gkr 2 5 



S—r 



> -^fcV + £iogfc. 

Tlius we have reduced the problem of bounding 1/Gr from below to solving a 
quadratic inequality; and then optimizing it for k. In the case 

(10.8) ^Xek^r^\ogk<l, 

we get by applying + x > 1 + x/A for < a; < 1: 



1 -1 + Jl + Xek^r"^ logfc/2 1 , , 

g; ^ — AP.V2 ^ 

In the case r^Xe < 1/2 wc choose a positive integer fc so that 



1 \V3 / X 1/3 

' < fc < 2 ' 



2r2A(ry " " \2r^Xe 

As fc^logfc < fc3/2, ([To:8l) holds. Inserting logfc > - log(2Aer2)/3 into ([TOJll 
completes the proof of Lemma [6l 
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